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test  of  Hq :  F  is  exponential,  versus  H^:  F  is  UBU,  but  not  exponential.  Our  test  is  based 
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suitably  normalized,  is  established.  Then  using  a  consistent  estimator  of  the  null  standard 
deviation  of  n1/2J^,  an  asymptotically  distribution-free  test  is  obtained.  Finally,  using 
tests  for  the  censored  and  uncensored  models  we  develop  a  measure  of  the  efficiency  loss 
due  to  the  presence  of  censoring. 
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Testing  Uhether  Mew  Is  Better  Than  Used 
With  Randomly  Censored  Data 

by 

Yuan-Yan  Chen,  ttyles  Hollander,  and  Naftall  A.  Lartgberg 

Abstract. 

A  life  distribution  F  Is  new  better  than  used  (fBU)  if  F(x  ♦  y)  <. F(x)F(y) 
for  all  x,  y  >.  0  (F  =  1  -  F).  Using  a  randomly  censored  sample  of  size  n 
from  F,  we  propose  a  test  of  HQ:  F  is  exponential,  versus  H.:  F  is  NBU, 
but  not  exponential.  Our  test  is  based  on  the  statistic 

3  //Fn(x  +  y)dFn(x)dFn(y),  where  Fn  Is  the  product  limit  estimator  of  F, 
Introduced  by  Kaplan  and  Meier  (1958). 

Under  mild  regularity  on  the  amount  of  censoring,  the  asymptotic  normality 
of  Jp,  suitably  normalized,  is  established.  Then  using  a  consistent  esti¬ 
mator  of  the  null  standard  deviation  of  n1^2,)^,  an  asymptotically  distribu¬ 
tion-free  test  is  obtained.  Finally,  using  tests  for  the  censored  and  un¬ 
censored  models  we  develop  a  measure  of  the  efficiency  loss  due  to  the  presence 
of  censoring. 

Key  words:  New  better  than  used,  exponential ity,  hypothesis  test, 
censored  data. 
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1.  Introduction  and  Sumary.  A  life  distribution  F  (a  distribution  function  .  ,  j 

r . . „  tL  t-lv**-***"  ' 

(d.f.)  such  that  F(x)  »  0  for  x  <  0),  with  survival  function's  1  -  ft,  qp 
Is  new  better  than  used  (HBUl.4f  }-P 

(1.1)  T{x  ♦  y)  <Jk1f(y.)  for  x,  y  <  [0,  •). 


S>The  dual  notion  of  a  new  worse  than  used  (NWU)  life  d.f.  is  defined  by  revers¬ 
ing  the  Inequality  .In-  (l-»4->£>  The  boundary  members  of  the  NBU  and  NUU  classes, 
obtained  by  Insisting  an  equal Ity^in  (lrHJ^are  the  exponential  d.f.'s. 

The  NBU  class  of  life  distributions  has  proved  to  be  very  useful  in 
performing  analyses  of  llfelengths.  These  d.f.'s  provide  readily  interpretable 
models  for  describing  wearout,  play  a  fundamental  role  in  studies  of  replace¬ 
ment  policies  (Marshall,  Proschan,  1972),  and  shock  models  (Esary,  Marshall, 
and  Proschan,  1973),  and  have  desirable  closure  properties  (c.f.  Barlow 
and  Proschan,  1975).^____ 

Hollander  and  Proschan  (HP)  (1972)  develop  a  test  of 


(1.2)  Hq *  F(x)  *  1  -  exp(-x/y),  x  >  0,  u  >  0  (u  unspecified), 
versus 

(1.3)  F  Is  NBU,  but  not  exponential. 


based  on  a  random  sample  X^,  ...,  Xn  from  a  continuous  life  distribution 
F.  The  hypothesis  Hq  asserts  that  a  new  Item  has  stochastically  the  same 
llfelength  as  a  used  Item  of  any  age,  where  the  alternative  H^  states  that 
a  new  Item  has  stochastically  greater  llfelength  than  a  used  item  of  any 
age.  The  HP  (1972)  test  Is  motivated  by  considering  the  parameter 


r(F)  *  /  f[F(x)V(y)  -  V(x  +  y)]dF(x)dF(y)  * 

(1.4)  0  0 

•  •  def 

*  1/4  -  /  / F(x  +  y)dF(x)dF(y)  =  1/4  -  a(F). 

0  0 

Viewing  y(F)  as  a  measure  of  the  deviation  of  F  from  exponential Ity  towards 
NBU  [or  NWU]  alternatives,  HP  (1972)  replace  F  by  Gn»  the  empirical  d.f. 
of  Xj,  ...»  Xn,  and  suggest  rejecting  Hq  in  favor  of  Hj  if  a(G0)  is  too 
small  [Hq  Is  rejected,  in  favor  of  H|t  F  is  NUU,  but  not  exponential,  if 
A(Gn)  Is  too  large.]  For  furtiier  details  about  the  test  see  HP  (1972), 

Hollander  and  VJolfe  (1973),  Cox  and  Hinkley  (1974),  and  Randles  and  Wolfe  (1979). 

In  this  paper  we  consider  a  randomly  censored  model  where  we  do  not 
get  to  observe  a  complete  sample  of  X's.  Let  Xj,  Xg,  ...  be  Independent  Identically 
distributed  (l.i.d.)  random  variables  (r.v.'s)  having  a  common  continuous 
life  d.f.  F.  The  X* s  represent  lifelengths  of  identical  items.  Let 
Yg,  •••  be  t*i*d.  r.v.'s  having  a  common  continuous  d.f.  H.  The  Y's 
represent  the  random  times  to  right-censorship.  Throughout  we  assune  the 
X's  and  Y's  are  mutually  Independent  and  the  pairs  (Xj,  Yj ),  (Xg,  Yg),  ... 
are  defined  on  a  common  probability  space  (a,  B,  P).  Further,  let  1(A) 
denote  the  indicator  function  of  the  set  A,  and  for  i  *  1,  ...,  n, 
let  Zj  ■  m1n{Xi$  Y^),  and  6^  -  I(Xj  <  Y^ ).  Based  on  the  incomplete  data 
set  (Zj,  Sj),  ...,  (Zn,  «n)  we  test  Hq,  given  by  (1.2),  against  Hj, 
given  by  (1.3).  The  censoring  d.f.  H  is  assumed  to  be  unknown  and  Is 
treated  as  a  nuisance  parameter.  Due  to  the  censoring,  tiie  empirical 
d.f.  GR  corresponding  to  F  cannot  be  computed.  Thus,  v/e  propose  to 
reject  Hq  In  favor  of  H^  for  small  values  of 

def  .  •  • 

“F")  ’  <  ’  l  fa'*  +  »*.<**.<»>• 


(1.5) 
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where  FR  is  the  Product  Limit  Estimator  (PLE)  of  F,  introduced  by  Kaplan 
and  Meier  (1958): 

_  def  .  «/n 

(1.6)  7(x>  -  1  -  F  (x)  -  JI  [(>>  -  1)(»  -  1  ♦  I)*1}  (,). 

"  "  0:  Z(1)<x) 

where  <  ...  <  Z^  denote  the  ordered  Z's,  Z^Qj  *  0,  «  1,  and 

&(l)»  •••»  6(n)  are  the  5's  corresponding  to  Zjjj,  ...»  Z^nj  respectively.  In 

(1.6) ,  we  treat  Z^nj  as  a  death  (whether  or  not  it  actually  Is)  so  that  *  1. 

Furthermore,  although  our  assumptions  preclude  the  possibility  of  ties.  In 
practice  ties  will  occur.  When  censored  observations  are  tied  with  uncensored 
observations,  our  convention,  when  forming  the  list  of  the  ordered  Z’s,  is  to 
treat  uncensored  members  of  the  tie  as  preceding  the  censored  members  of  the  tie. 

For  computational  purposes,  it  is  convenient  to  write  £  as: 


-  I?.iI7.itC  “  {(n  -  k)(n  -  fc  ♦  l)'1)  (lt,:l 

"  ,lJ1  *«  z(k)iz(l)«U)» 

1 ”1  1  ilr\  J-l  I  ii.l 

C  n  {(n  -  r)(n  -  r  +  lp)  'r'][  n  t(n  -  s)(n  -  s  +  l)"1)  (s}] 


5=1 


(n  -  i  +  l)”1(n  -  j  +  1)’lfi(i)6(j)K 

In  clinical  trials  the  X's  may  be  times,  measured  from  date  of  diagnosis 
of  a  disease,  to  relapse.  NBU  alternatives  may  be  preferable  to  increasing 
failure  rate  alternatives  because  the  latter  class  insists  upon  a  non-decreas¬ 
ing  failure  rate  whereas,  in  this  medical  context,  we  may  expect  the 
failure  rate  to  increase  (at  least  for  a  short  period  of  time)  after 
treatment  begins.  Incomplete  observations  can  arise  at  the  time  of  data 
analysis  due  to,  for  example,  dropout  or  patients  who  have  not  yet 
relapsed.  In  this  situation,  it  is  appropriate  to  use  to  test  Hq  vs.  Hj. 

Marshall  and  Proschan  (1972)  consider  age  replacement  policies  and 
block  replacement  policies.  Under  an  age  replacement  policy,  a  unit  is 
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replaced  upon  failure  or  upon  reaching  a  specified  age  T,  whichever  comes 
first.  Under  a  block  replacement  policy,  a  replacement  Is  made  whenever  a 
failure  occurs,  and  additionally  at  specified  times  T,  2T,  3T,  ...,  . 

Marshall  and  Proschan  show  that  a  necessary  and  sufficient  condition  for  failure- 
free  Intervals  to  be  stochastically  larger  (smaller)  under  age  replacement  than 
under  a  policy  of  replacement  at  failure  only  Is  that  the  underlying  distribution 
be  NBU  (NWU).  Marshall  and  Proschan  also  show  that  a  necessary  and  sufficient 
condition  that  the  nunber  of  failures  in  a  specified  CO,  t]  be  stochastically 
smaller  (larger)  under  age  replacement  than  under  a  policy  of  replacement  at 
failure  only  is  that  the  underlying  distribution  be  NBU  (fWU).  Similar  compari¬ 
sons  hold  for  block  replacement.  Thus  in  reaching  a  decision  as  to  whether  to 
use  an  age  (block)  replacement  policy  or  not,  it  is  important  to  investigate 
whether  or  not  the  unaerlylng  distribution  Is  NBU.  If  lifelength  times  are 
censored,  the  test  based  on  facilitates  such  an  investigation. 

Other  references  describing  situations  where  it  is  important  to  know 
whether  the  underly  distribution  Is  1CU  are  Esary,  Marshall  and  Proschan 
(1973)  In  the  context  of  shock  models,  and  El-Neweihi,  Proschan,  and  Sethuraman 
(1978)  In  the  context  of  multiple  coherent  systems. 

In  Section  2  we  establish  the  asymptotic  normality  of  the  sequence 

n^2{,£  -  a(F)}  under  the  assumptions: 
n 

(A.l)  The  supports  of  F  and  H  are  equal  to  [0,  •), 

(A.2)  sup{(T(x)]1"€Cff(x)]'1,X€[0,-)}  <  - 

for  some  nonnegative  real  number  c. 
and 

(A.3)  The  processes  {n^2(Fn(t)  -  F(t)),  tc(-»,»)>  converge  weakly 
to  a  Gaussian  process  with  mean  zero  and  covariance  kernel  given  by  (2.1). 
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Condition  (A.2)  restricts  the  amount  of  censoring  allowed  in  the  model. 

To  see  this  In  a  simple  case,  consider  the  proportional  hazards  model  where 
W  »  fF)6  for  some  3  >  0.  Then  PCX^  Y^J  *  (3  +  1)“*,  and  condition  (A.2) 

implies  that  3  £  1.  Thus,  in  the  proportional  hazards  model,  the  Jjj  test 

Is  Inappropriate  when  the  expected  amount  of  censoring  P(Y^  <  Xj)  exceeds 
50*. 

The  null  asymptotic  mean  of  is  1/4,  Independent  of  the  nuisance 
parameters  u  and  H.  However,  the  null  asymptotic  variance  of  n^2J?  does 

VI 

depend  on  v  and  H  and  must  be  estimated  from  the  data.  A  consistent  estl- 

mator,  a|j,  given  by  (3.3),  is  derived  in  Section  3.  The  approximate  o-level 

test  rejects  Hq  Is  favor  of  if  n^2{J^  -  (1/4  ^o'1  <  -  z^,  where  zo  is 
the  upper  a-percentile  of  a  standard  normal  distribution.  In  Section  3 
we  also  show  that  this  asymptotically  distribution-free  test  is,  under 
suitable  regularity,  consistent  against  all  continuous  NBU  alternatives. 

Section  4  develops  a  measure  of  the  loss  in  efficiency  due  to  the  presence 
of  censoring.  This  measure  Is  derived  using  the  HP  (1972)  ?»U  test  and 
its  generalization  herein  proposed  based  on  J^.  In  certain  instances,  this 
measure  assumes  values  close  to  P(Xj  <  Yj);  the  latter  also  being  a  (rough) 
measure  of  the  loss  of  Information  due  to  censoring. 

Section  5  contains  an  application  of  the  statistic  to  some  survival 

data. 

2.  Asymptotic  Normality  of  the  NBU  Test  Statistic.  In  this  section  we 

establish  the  asymptotic  normality  of  the  test  statistic  «£,  defined  by  (1.5). 

n 

Let  F(t)  -  F(t)lT(t),  tc(— ,-),  and  let  (4(t),  te(-»,«)>  be  a 
Gaussian  process  with  mean  zero  and  covariance  kernel  given  by: 


^(tjrtsj/c^zjHzjr^z),  o  <  s  <  t<  -, 

(2.1)  E*(t)*(s)-^  0 

{0  ,  s  <  0  or  t  <  0. 

Unless  otherwise  specified,  all  limits  are  evaluated  as  n  ♦  •,  and  all  Inte¬ 
grals  range  over  (-•,  -). 

First  we  state  the  main  result  of  this  section. 

Theorem  2.1.  Assume  that  conditions  (A.l),  (A.2),  and  (A.3),  given 
in  Section  1,  hold.  Then  n^2{J^  -  a( F )>  converges  in  distribution  to  a 

p 

normal  r.v.  with  mean  zero  and  variance  a  ,  given  by: 

(2.2)  o2  -  fff/ElMt  +  s)  -  2*(t  -  s)][*(u  +  v)  -  2*(u  -  v)])dF(s)dF(t)dF(u)dF(v). 
Note  that  for  n  *  1,  2,  ...» 

"  a(F)  -  </JtFn(x  +  y)  -  F(x  ♦  y)]dFn(x)dFn(y) 

-  //CFn(x  +  y)  -  F(x  +  y)]dFn(x)dF(y)} 

+  t//CTn(x  +  y)  -  F(x  +  y)]dFn(x)dF(y) 

-  //[F^x  +  y)  -  F(x  +  y)]dF(x)dF(y)} 

+  //!Tn(x  +  y)  -  F(x  +  y)]dF(x)dF(y) 

+  //F(x  +  y)dFn(x)dFn(y)  -  // F(x  +  y)dF(x)dFn(y) 

♦  //F(x  +  y)dF(x)dFn(y)  -  ffF(x  *  y)dF(x)dF(y). 

Upon  integration  by  parts  and  change  of  variable  we  obtain  that: 


//F(x  +  y)dFn(x)dFn(y)  -  //F(x  +  y)dF(x)dFn(y) 

■  -  //CFfl(x  -  y)  -  F(x  -  y)]dFn(x)dF(y),  n  *  1,  2,  .... 


-  /  - 


and  that 


f/T(x  +  y)dF(x)dFn(y)  -  JJF(x  +  y)dF(x)dF(y) 

=  -//lTn(x  -  y)  -  F(x  -  y )]dF(x)dF(y),  n  =  1,  2 . 


Thus  for  n  *  1,  2,  ... , 


n1/2{JCn  -  A(F)}  =  Bn  .  +  Bn  9  -  Bn  ,  +  Bn  ,, 
n  n,i  n,2  n»3  n,4* 


where 


Bn,l  "  //nl/2C'rn(x  +  y)  ~  F(x  +  y)]dFn(x)dFn(y)  - 

-  //n1/2[Fn(x  +  y)  -  F(x  +  y)]dFn(x)dF(y), 

Bn,2 3  i7nl/ziTn(x  +  y>  -  ^x  +  y)]dFn(x)dfr(y)  - 

-  //-1/2CFn(x  +  y)  -  F(x  +  y)]dF(x)dF(y), 

Bn,3  “  //nl/2C^(x  -  y)  -  F(x  -  y)]dFn(x)dF(y)  - 

-  J/n1/2CFn(x  -  y)  -  F(x  -  y)]dF(x)dF(y), 

and 


Bn,4  =  //"1/2Fn(x  +  y)  -  F(x  +y)  -  2{FR(x  -  y)  -  F(x  -  y)}]dF(x)dF(y). 


Consequently  to  prove  the  result  of  Theorem  2.1  it  suffices,  by  Slutsky's 

Theorem  [Billingsley  (1968),  p.  49],  to  show  that  B„  .»  B_  ,,  and  B„  ,  converge 

n»i  n  »c  n 

in  probability  to  zero,  and  that  B  .  converges  in  distribution  to  a  normal 

n 

r.v.  with  mean  zero  and  variance  o  ,  given  by  (2.2). 

First  we  prove  that  Bn  B_  ,,  and  B„  ,  converge  in  probability  to 

f)f  i  n  1 1  n 

zero,  and  that  ^  converges  In  distribution.  Then  we  prove  that  the  limit- 

*  2 
ing  d.f.  of  B_  .  Is  normal  with  mean  zero  and  variance  o  . 

To  establish  the  convergence  of  Bn  j  through  Bn  ^  we  Introduce  a  notation 

and  four  lemmas.  Let  0  ®  (<j»:  <»  Is  real  valued,  bounded,  and  right-continuous 
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function  defined  on  (-•,  •),  with  finite  left-hand  limits  at  each  t«(— ,»), 
and  finite  limits  at  t  «  ±  •}.  Throughout  we  view  D  as  a  metric  space  with 
the  Skorohod  metric  [Billingsley  (1968),  p.  112]. 

lemma  2.2.  Let  <  0,  and  ye(-«,«).  Then  f\ |»(x  +  y)dFn(x),  and 
/4»(x  -  y)dFn(x)  converge  w.p.l,  to  f\ |»(x  +  y)dF(x),  and  to  /< |>(x  -  y)dF(x), 
respectively. 

Proof.  There  is  a  set  e  8,  =  1,  such  that  Fn(x,  <u)  converges 

to  F(x)  for  every  x«(«,«),  and  <a  c  [c.f.  Peterson  (1977),  Th.  3.3, 
or  Langberg,  Proschan,  and  Quinzi  (1980),  Th.  4.9].  Note  that  the  sets 
of  discontinuities  of  the  functions  *>(•  +  y),  and  4»( *  -  y)  are  countable 
[Billingsley  (1968),  p.  110].  Since  F  is  continuous  these  sets  have  F-measure 
zero.  Consequently  thi  desired  results  follow  by  the  Helly-Bray  Lemma  [Breiman 
(1968),  p.  163,  Th.  8.12].  (| 

Lemma  2.3.  Let  ^  «  0.  Then  Jj>(x  +  y )dFn(x)dF(y),  and  JjXx  -  y )dFn(x)dF(y), 
converge  w.p.l  to  JJ<|i(x  +  y )dF(x)dF(y),  and  ff\ j»( x  -  y )dF(x)dF(y),  respectively. 

Proof.  Note  that  /i|»(x  +  *)dFn(x),  and  /i»(x  -  *)dFn(x)  are  sequences 
of  bounded  functions.  By  Lemma  2.2  these  sequences  converge  w.p.l  to 
/i|»(x  +  «)dF(x),  and  f\ |>(x  -  •  )dF(x), respectively.  Consequently  the  desired 
results  follow  by  the  Dominated  Convergence  Theorem.  || 

Lemma  2.4.  Let  <|»  be  a  continuous  function  in  D.  Then 
//♦(x  +  y)dFn(x)dFn(y)  -  /J>(x  +  y)dFn(x)dF(y),  and  //t(x  -  y )dFn(x)dFn(y)  - 
//♦(x  -  y)dFn(x)dF(y),  converge  w.p.l  to  zero. 

Proof.  To  prove  the  desired  results  it  suffices,  by  Lemma  2.3,  to 
show  that  fU(x  +  y)dFn(x)dFn(y),  and  //*(x  -  y)dFn(x)dFn(y),  converge  w.p.l- 
to  //i»(x  +  y)dF(x)dF(y),  and  to  ff\ |/(x  -  y )dF(x)dF(y),  respectively. 


We  now  prove  the  preceding  t wo  statements.  There  is  a  set  «  rj , 

P{ftj}  *  1,  such  that  Fn(x,  u)  converges  to  F(x)  for  all  X6(-»,«),  and 

a)  c  flj.  Consequently  the  ^recedin  ',  two  statements  follow  by  the  Helly-Bray 

Lenina  [Billingsley  (1968),  p.  11,  Th.  2.1  (11)].  || 

Lemma  2.5.  Assume  (A.l)  holds.  Then  the  Gaussian  process 

4(t),  t€(~»,«)},  with  mean  zero  and  covariance  kernel  given  by  (2.1), 

has  continuous  paths  w.p.l. 

t  , 

Proof.  Let  g(t)  =  jlK(z)F(z)]  1dF(z),  tcCO,-).  By  (A.l)  and  the  continuity 
0 

of  F,  g  is  strictly  increasing  and  continuous.  Let  $j(t)  be  a  stochastic 
process  given  by: 


^(t) 


fcW^tJxrMg’^t)).  te[o,«), 

J 

I  o  ,  te(-»,0). 


Clearly  {^(t),  tct-,-)}  is  a  Gaussian  process  with  mean  zero,  and 
covariance  kernel  given  by: 


E*1(t)^1(s)  = 


0  <  s  <  t  <  ® 
s  <  0  or  t  <  0. 


Thus  Uj(t),  tc[0 »• ))  is  a  standard  Wiener  process. 

Note  that  g"*  is  continuous  and  strictly  increasing  to  «,  g_1(0)*  0, 
and  that  {^(t),  t«(-»,«)}  has  continuous  paths  w.p.l  [Breiman  (1968), 
p.  257],  Consequently  the  desired  result  follows  from  the  definition  of 
(<frj(t),  t«(— ,-)).  || 

Breslow  and  Crowley  (1974)  and  Peterson  (1977)  prove  that  the  processes 
{n*^{Fn(t)  -  F(t)>,  t€(— ,T)}  converge  weakly  to  the  process  ($(t),  t«(-»,  T )) 
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for  all  Te( -•.•),  provided  (A.l)  holds.  To  prove  that  Bn  Bn  and  B„  , 

II  U  jC  llfv 

converge  in  probability  to  zero,  and  that  Bn  ^  converges  in  distribution 

*  o 

to  a  normal  r.v.  with  mean  zero  and  variance  a ,  we  must  assume  that  Breslow 

and  Crowley's  result  holds  for  T  =  «.  Since  (A.2)  restricts  the  amount 

of  censoring  allowed  in  the  model  we  conjecture  that  under  (A.l),  (A.2), 

(n^2{Fn(t)  -  F(t)},  t£(— ,»))  converges  weakly  to  ($(t),  tc(-»,»)}.  Me 

assume  this  conjecture  to  prove  Lemmas  2.6,  2.7  and  2.9  that  follow. 

Now  we  proceed  to  prove  that  Bn  B_  ,,  and  Bn  ,  converge  in  probability 

to  zero,  and  that  Bn  ^  converges  in  distribution.  A  lemma  is  needed. 

Lemma  2.6.  Assume  (A.3)  holds.  Then  the  processes  {n^2{F_(x  +  v)  - 
"  n 

F(x  +  y)  -  2[F^(x  -  y)  -  F(x  -  y)]>,  x,  ye(— converge  weakly  to  the 

process  t$(x  +  y)  -  2y(x  -  y),  x,  ye(— ,«)). 

2 

Proof.  Let  D  *  {<  >,  i|»2  e  Dl,  be  a  mertic  space  with  the 

metric  induced  by  the  one  of  D.  By  a  standard  argument  the  bivariate  processes 
(n1^2  <  Fn(t)  -  F(t),  Fn(s)  -  F(s)  >,  t,  se(— ,»)>  converge  weakly  to  the 
bivariate  process  {<  $(t),  $(s)  >,  t,  s€(-»,«)}.  Thus,  by  the  Continuous 
Mapping  Theorem  (Billingsley  (1968),  P.  30,  Th.  5.1)  the  processes 
{n^2{Fn(t)  -  F(t)  -  2[Fn(s)  -  F(s)]>,  t,  S€(-»,»)>  converge  weakly  to  the 
process  U(t)  -  2$(s),  t,  se(-»,-)}.  Consequently  the  desired  result  follows.  || 

Ue  no  vjstaylis.j  The  convergence  of  Bp  ^  through  Bn  Some  notation 

12  12  * 
is  useful.  Let  Q  ,  Q  ,  Q*,  Qj|,  be  the  probability  measures  on  D  induced  by 

the  processes  4(t),  t«(— ,-)},  {<fr(x  +  y)  -  2$(x  -  y),  x,  ye(— ,•)), 
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{n^2{7n(t:)  -  F(t)},  and  (n1/2<Fn(x  ♦  y)  -  F(x  +  y)  - 

2lTn(x  -  y)  -  F(x  -  y )]> ,  x,  ye(— ,  «)},  n  =  1,  2,  ....respectively.  Let 
Sj,  Sg  be  two  sets,  let  A  be  a  subset  of  S2  and  let  £  be  a  mapping  from  Sj 
to  S2;  then  e_1(A)  *  {s:  s  e  S^,  5(s)  e  A). 

Lemma  2.7.  Assume  (A.l),  (A. 2),  and  (A. 3)  hold.  Then: 

(a)  Bn  Bn2,  and  Bn  3  converge  in  probability  to  zero, 
and 

(b)  B  .  converges  in  distribution  to  the  r.v.  //[♦(x  ♦  y)  - 

n 

2<fr(x  -  y)]dF(x)dF(y). 


Proof.  For  e  D,  and  n  =  1,  2,  ...  let 


$n>1(*)  a  JJ<Hx  +  y)dFn(x)dFn(y)  -  //*(x  +  y)dFn(x)dF(y), 

€„,*(♦)  »  //♦(x  +  y)dFn(x)dF(y)  -  //^(x  +  y)dF(x)dF(y), 

Ln>3U)  »  //.(x  -  y)dFn(x)dF(y)  -  JJ<Kx  -  y)dF(x)dF(y), 
and 

C4(t|>)  =  ffMx  +  y)  -  2xp(x  -  y)]dF(x)dF(y). 

1  1  2 
The  probabilities  Q*  converge  weakly  to  Q  .  By  Lemma  2.6  Qn  converges 

weakly  to  Q2.  By  Lemma  2, St  the  supports  of  Q*  and  Q2  coincide  with  the 

set  of  all  continuous  functions  in  D.  By  the  definitions  of  the  mappings 

and  probability  measures: 

qU-JM-,  x3>  =  PCB  .  <  x},  xc(— q  =  1,  2,  3,  n,  *  1,  2,  .... 

n  n,q  n ** 

Q^C”1  {(— .x])  *  <  x),  xc(— ,»),  n  =  1,  2,  ... 

and 

qV1^— ,  x]J  ■  P(/J[$(u  +  v)  -  2${u  -  v)]dF(u)dF(v)  <  x),  x«(-«,-). 
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Thus  to  obtain  the  desired  results  it  suffices  to  show,  by  the  Extended 
Continuous  Mapping  Theorem  [Billingsley  (1968),  p.  34,  Th.  5.5],  that  for 
every  sequence  y  0  that  converges  to  a  continuous  function  f  «  D, 
lim  5n>q(*n)  *  °*  "’P*1  for  q  “  1§  2»  3»  and  1im  *^n^  a 

Me  o  pro v:  the  preceding  statements.  Let  y  D,  n  •  1,  2 . 

ud  let  f  be  a  continuous  function  in  0.  Assume  lim  =  ♦.  By  a  well-known 
result  [Billingsley  (1968),  p.  112]: 

lim  sup{|*n(x)  -  <J»(x){,  x€(-»,»))  «  0. 

By  Lemma  2.4  lim  £  .(*)  »  0,  w.p.l.  By  Lemma  2.3  lim  £_  „(*)  *  0,  q  *  2,  3. 

Consequently  by  simple  integral  evaluations  we  obtain  that  lim  tn  q(fn)  *  0. 

w.p.l.  for  q  *  1,  2,  ?,  and  that  lim  6(*n)  *  ?(♦).  )| 

Vie  are  ready  to  show  that  the  limiting  d.f.  of  Bn  ^  is  normal  with 

mean  zero  and  variance  a2,  given  by  (2.2).  First,  we  show  that  under  (A.2), 

2  . 

0  V  •# 

Lemma  2.8.  Assume  (A.2)  holds.  Then  a2,  given  by  (2.2),  is  finite. 

Proof.  Note  that  for  a,  b€(-»,»),  (a  -  b)2  <  2(a2  +  b2).  Thus,  by 
the  Cauchy-Schwartz  Inequality: 

o2  <  //£[♦( s  +  t)  -  2f(s  -  t)]2dF(t)dF(s)  < 

<  2//[E4(s  +  t)}2  +  4E(*(s  -  t))2]dF(s)dF(t) 

<  10  sup{E{$(t)}2,  t«[0,-)}. 

Hence  to  prove  the  desired  result  it  suffices  to  show  that 

supflTtt)]2/1 CK(z)F(z)]-1dF{i),  By  (A.zj.tltji)]'1  <  cCTti)]**1 

/  0 

for  all  z«[0,«),  some  c«{0,*»)»  and  some  nonnegative  real  number  «.  Consequently 
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dWAkUFU)]-1^*)  icCRt^/Vutf^dFU).  tc[0.-). 

0  0 

Since  F  Is  continuous 

t  \U  -  2)_1{CF(t)]2  -  CF(t)]*>.  e  *  2. 

CT(t)]2/  CF<z)]*’3dF(z)  =  (  for  t«[0,-  . 

0  j£F(t)]2  -  U(t)]2tnF(t)  ,  c  -  2, 

The  desired  result  follows  now  by  simple  limiting  arguments.  || 

Finally,  we  Identify  the  limiting  d.f.  of  Bn 

Lemma  2.9.  Assune  (A.1),(A.2),  n-1  (A.’)  hold.  T  'i  ccnvonos  In  distri- 
button  to  a  normal  r.v.  with  mean  0  and  variance  a ,  given  by  (2.2). 

Proof.  By  Lemma  2.7  (b)  it  suffices  to  show  that  the  r.v.  JJ[.(x  +  y)  - 

O 

24 (x  -  y)]dF(x)dF(y)  Is  normal  with  mean  zero  and  variance  o  .  since  the 
process  4(x  +  y)  -2*(x  -  y),  x,  ye(-»,»)}  Is  Gaussian  the  desired  result 
follows  by  the  theory  of  stochastic  Integration  [c.f.  Parzen  (1962),  p.  78].  || 

3.  Consistency.  Let  f(z)  *  z3[l  +  4*nz  +  4(*nz)2]/16,  0  <  z  <  1,  and  *  0 
for  z  *  0,  let  u  *  /xdF(x),  and  let  n  *  [P(Xj  <  Further,  let 

Up  *  ^(*1^1*  Kp(t)  =  n  Ej.jlCZj  >  t),  n  *  1,  2,  ...»  tc(-»,»). 

Finally,  let 

(3.1)  o2(e)  *  JfUJCfcHtnz)]"^,  B«(0,«), 

0 


(3.3)  %  »  S2^),  n  «  1,  2 . 

2 

For  computational  purposes  dp  can  be  written  as 


S2  -  (128 1*1  ♦  t";Jn(n  -  1  t  ir'(n  -  i)"1C(128)‘1  - 

(3ar1z(t)«nr1  *  (i6r1z^)(inr-j«pMz(t)(i:„r1)]  - 
ntazar1  -(32)',z(n)(;n)-1  +  (i6)'lz^;n)'2: 
exp{-4Z(rl)(5nr‘>- 

In  this  section  we  show  that,  under  HQ,  S2  Is  a  consistent  estimator 
of  a2  provided  o^(0)  Is  finite  In  a  neighborhood  of  y.  He  then  show,  under 
(A.l),  (A.2)  and  (A.3),  and  the  assunptions:  y  <  »  and  <r2(o >  <  •  in  a  neigh¬ 
borhood  of  n,  that  the  approximate  a-level  test,  which  rejects  Hq  In  favor 
of  Hj  If  n1,2{jjj  -  (1/4  JJo"1  <  -  za,  is  consistent  against  all  continuous 
MBU  alternatives.  He  conclude  the  section  by  presenting  a  sufficient  condition 
for  o2(e)  to  be  finite  at  M0,«). 

Mow  we  show  the  consistency  of  ojj  under  Hg.  The  proof  of  consistency 
uses  several  lemmas.  We  first  show  that  ar2(e )  converges  in  probability 
to  o2(e),  provided  o2(e)  <  •.  Then  we  show  that  o2  converges  in  probability 
to  o2(n),  provided  o2(e)  is  finite  in  an  interval  containing  n.  Finally, 
using  the  previous  results,  we  obtain  the  consistency  of  under  Hg. 

To  show  that  o2(G)  converges  to  o2(e)  we  need  a  well-known  proposition, 
stated  for  the  sake  of  completeness,  and  a  lemma. 

Proposition  3.1.  [David  (1970),  p.  18]  Let  U^j  <  ...  <  U^nj  be  the 
order  statistics  of  a  sample  of  size  n,  n  *  2,  3,  ...,  taken  from  a  continuous 
d.f .  6.  Further  let  u«(1nf{s:  G(s)  >  0),  •),  and  Gu(t)  -  GCtJCGfu)]"1, 
t«(-»,  ).  Then  the  conditional  random  vector  {(U^j,  •  ••>  u(n-l)^u(n)  " 
is  stochastically  equal  to  the  order  statistics  of  a  sample  of  size  n  -  1 
taken  from  the  d.f .  s  . 
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Laima  3.2.  Let  MO,**),  and  MO,l).  Then 

E/f(z)Cln(-8inz)]“1I(-tnz  <  8“lZ^nj)dz  < 

6  , 

<  3/f(z)[T<(-8inz)]  *dz,  n  *  1,  2,  ...  . 

0 

Proof.  Note  that  for  n  =  1,  2,  ..., 

/ 

/ 

«  i  i 

E/f(z)[Kn(-e*nz)]  Mf-tnz  <  e  ‘Z^Jdz  = 

*  /f(zH/!einzE{C1Si(-em)]  lZ(n)  *  «WPtt(n)  <  u>Mz. 

Now  let  zc(0,l),  u«(-etnz,»),  q(u,  z)  *  [K(u)]_1K(-0)tnz),  and 
p(u,  z)  *  1  -  q(u,  z).  By  Proposition  3.1  the  conditional  r.v. 
{[^(-einzjr1)!^^  *  u}  Is  stochastically  equal  to  n[B(n  -  1,  p(u,  z))  +  lj"1, 
where  B(n  -  1,  p(u,  z))  Is  a  binomial  r.v.  with  parameters  n  -  1,  and 
p(u,  z).  Thus 

E{C^(-8tnz)]”1jZ^nj  »  u}  » 

■  n  ♦  l)'l(J'l)[p(».  z)]J[q(u, 

•  CpOi,  I)r'(l  -  [q(u.  z)Jn). 

Let  m(n)  -  n  for  n  *  2,  4,  6,  ....  and  »  n  +  1  for  n  *  1,  3,  5 . 

Then 


cp(u.  ,)r*<i  -  cq<».  i)d")  <  cp(».  zir'u  -  c«cu.  ^)^*(",) 

<  2[p(u,  z))*1!!  -  [q(u,  z)f<">/2)  -  »$>«[,<».  r)]\  n  -  l,  2 . 
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By  direct  evaluation  we  obtain  that: 

7  E{CKn(-»*n2)r1|Z(  ,  *  u)dP(Z,  v  <  ul 

-etnz  "  v  '  17 

<  icflj^n  -  <  3[I(-e*nz)]  ’l . 

Consequently  the  desired  result  follows.  || 

o  O 

We  are  ready  to  prove  that  3^j(e)  converges  In  probability  to  <s  (e). 

Lemma  3.3.  Let  6e(0,<»).  Assume  o2(0),  given  in  (3.1),  Is  finite.  Then 

p  -  11m  o2(e)  *  o2(0). 

Proof.  Let  X,  6«(0,»).  Then 

P{|o2(e)  -  o2(e)|  >  x>  <  P(/f ( z )[Kn( -0inz )]“* I ( - tnz  <  0_1Z(n))dz  >  x/3) 

«  . 

+  Pt/ffzJCKf-etnzJl'^dz  >  X/3) 

0 

+  P(/f ( z ) | [ Kn( -8tnz )]“*!( -tnz  <  -  C"K(-e*nz )3“1|dz  >  X/3). 

By  the  G1 i v enko- Cant ell 1  Lemma: 

1  .  •  1  .1 
Hm  /(fUXl^-MnzH^K-ftnz  <  «  Z(n))d2  -  /f(z)[K(-e*nz)]  Adz.  w.p.l. 

Thus  by  Lemma  3.2,  and  the  Chebyshev  Inequality: 

Urn  P{|d2(e)  -  o2(e)|  >  x)  <  6x"1/f(z)C?(-etnz)]“1dz. 
n  0 

Now  since  o2(0)  <  • 

«  , 

11m.  Jf(z)Ci<(-®inz)]“1dz  «  0. 

«-K)  0 
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Consequently  the  desired  result  follows.  || 

aO  O 

We  now  prove  that  converges  In  probability  to  o  (n). 

O 

Lenina  3.4.  Assume  that  n  <  •  and  that  a  (0)  is  finite  in  an  open 
interval  that  contains  n.  Then 

*2  2 

p  -  lim  o|j  *  o'  (n ). 

Proof.  Let  0  <  6  <  Sq  <  n,  Xc(0,«»),  and  let  AR  »  {«:  |nn(»)  -  n|  >6). 
Assume  o2(0)  <  •  for  e€[n  -  n  +  6Q3.  By  the  monotonicity  of  o2(0)  in 
0  for  n  *  1,  2 . we  obtain  by  some  simple  computations  that: 

Pt|o2  -  o2(n)|  >  X}  <  P{|o2(n  ♦  6)  -  ff2(n  +  «)|  >  X/4} 

+  P(|o2(n  -  6  )  -  o2(n  -  i)|  )  A/i)  +  4x  ^|o2(n  +  6)  -  o2(n)| 

+  4x’1Jo2(n  -  6)  -  <£(n)|  +  P{An),  n  *  1,  2.  ...  . 

By  the  Weak  Law  of  Large  Numbers,  1  im  poy  =  o. 

Thus  by  Lemma  3.3: 

Tim  P{|S2  -  o2(n)|  >  X)  <  4k’1|o2(u  ♦  6)  -  o2(n)!  +  4x*1|o2(n  -  6)  -  o2(n)|. 

Consequently  the  desired  result  follows  from  the  continuity  of  ^(e )  in 

[n  -  «0,  n  +  4q],  by  letting  6  ♦  0+,  || 

2 

We  obtain  now  the  consistency  of  o*  under  Hq. 

2 

Theorem  3.5.  Assume  o  (e)  is  finite  In  an  Interval  that  contains  u* 

Then  under  Hq 

o  o 
p  -  11m  St  ■  o  . 
n 

Proof.  Note  that  under  HQ,  »  ■  n  and  that  o2  ■  <?(n).  Consequently 
the  desired  result  follows  by  Lemma  3.4.  || 
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Next  v<e  show  that  our  test  is  consistent. 

Theorem  3.6.  Assume  (A.l),  (A.2),  and  (A.3)  hold.  Further,  assume 

that  u  <  »,  and  that  o2(0)  is  finite  in  an  interval  that  contains  n.  Then 

the  test,  which  rejects  Hn  in  favor  A  II.  if  n^2{J^  -  (1/4)>5*1<  -2  ,  is 

0  1  n  n  o 

consistent  against  all  continuous  NBL'  alternatives. 

Proof.  Note  that 

-  (1/4))^’  <  i(F)>  <  -zja  *  ill/2((l/4)  -  4(F))}, 

that  under  H,  (1/4)  •  a(f)  >  ).  an)  that  by  Lemna  3.4  p-  lim  am  *  o(n)  <  •. 
Consequently  the  desired  result  follows  by  Theorer,  2.1.  || 

Finally,  we  present  a  sufficient  condition  for  oZ(0)  <  •  at  0c(O,-). 

Lemma  3.1.  Let  ar(Q,<»).  resume 

(3.4)  lim  ^"^(-etnz)]"1  <  ®  for  some  pe(0,-0» 

z-*0 
o 

Then  o  (0)  <  «. 

Proof.  To  obtain  i:l:a  desired  result  is  sufficies  to  show  that  there 
is  a  4c(0,l),  such  chct 

c  _  1 

/f(z)[K(-0Jini;)?“idz  <  «. 

0 

We  show  no \i  the  preceding  inequality.  There  is  a  $«(0,1)  and  a  d«(0,»), 

such  that  z4-0  <OK(-eznz),  zclO.S).  Thus: 

<S  6 

jf{z )[](( -e « ru 4d  z  <  d/f( z)z0-4dz. 

C  0 

«  .  , 

Consequently  the  desired  result  follows  by  an  evaluation  of  /f(z)zB“^dz.  || 

0 


4.  Efficiency  loss  due  to  censoring.  Recall  that  the  test  is  a  generalization 
of  the  HP  (1972)  test  for  the  uncensored  model  based  on  the  statistic  Jn 
(see  equation  (1.5)  of  HP  (1972)).  In  this  section  we  study  the  efficiency 
loss  due  to  the  presence  of  censoring  by  comparing  the  power  of  the  test 
based  on  n  observations  in  the  uncensored  model  with  the  power  of  the  jjj 
test  based  on  n*  observations  in  the  randomly  censored  model. 

Let  F  be  a  parametric  family  within  the  NBU  class  with  F„  being  exponential 
Y  Yq 

with  scale  parameter  1  (for  example,  one  such  family  is  the  Ueibull 
F^(x)  ■  1  -  exp(-(x)Y),  y  >.  1  and  Yq  =  1 )  nv!  nssu -  o  t V.  rando  ly  censored  ioiiel  with 
F  *  F^  and  with  censoring  distribution  H.  Consider  the  sequence  of  alterna¬ 
tives  Yn  =  Yg  +  cn-1^,  with  c  >  0.  Let  8n(Yn)  be  the  power  of  the  approximate 
a-level  Jn  test  based  on  n  observations  in  the  uncensored  model  and  let 
Bn*(Yn)  denote  the  power  of  the  (approximate)  a-level  test  based  on  J®  for 
n*  observations  in  the  randomly  censored  model.  Consider  n*  «  h(n)  such 
that  11m  Bn(Yn)  31  lim  Bn*(Yn),  where  the  limiting  value  is  strictly  between 
Q,  and  1.  and  let  \  =  11m  n/n*.  The  value  of  k  can  be  viewed  as  a  measure 
of  the  efficiency  loss  due  to  censoring.  The  value  of  k  Is  adopted  from 
Pitman's  (cf.  Noether,  1955)  measure  of  asymptotic  relative  efficiency  but 
the  interpretation  of  k  must  be  modified  because  Jn  and  J®  are  not  competing 
tests  which  are  both  applicable  in  the  randomly  censored  model.  Roughly 
speaking,  for  large  n  and  NBU  alternatives  close  to  the  null  hypothesis 
of  exponential Ity,  the  test  requires  n/k  observations  from  the  randomly 
censored  model  to  do  as  well  as  the  Jn  test  applied  to  n  observations  from 

the  uncensored  model.  It  can  be  shown  that  since  JL  and  have  the  same 

n  n 

asymptotic  means,  k  reduces  to 


T 


■°r 
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def  c  9 

(4.1)  k  =  eH(Jc,  J)  «  (5/432 )/c*(l) 

1/2  2 

where  (5/432)  is  the  null  asymptotic  variance  of  n  Y  J  and  o  (1),  given 
by  (3.1),  is  the  null  asymptotic  variance  of  n^^J jj.  Thus  note  that  k  depends 
only  on  the  censoring  distribution  H,  and  not  on  the  parametric  family  Fv 
of  NBU  alternatives.  Hence  we  use  the  notation  e^(Jc,  J),  rather  than  ep  H(JC»  0: 
in  (4.1). 

lie  consider  the  cases  (i)  where  the  censoring  distribution  is  exponential, 

TTj(x)  *  1  for  x  <  0,  Hj(x)  *  exp(-Ax),  x  >  0,  and  (11)  where  the  censoring 

distribution  is  piecewise  exponential,  Tig(x)  =  1  for  x  <  0,  and  for  y  -  1 . . 

TT2(x)  -  cr  exp(-xrx),  s^  <  x  £  sr,  and  H2(x)  *  cm+1exp(-Vlx)»  sm  <  x  '*here 
r-1 

cp  *  exp  (  -I  Aj(Sj  -  -  1)  +  xrsr-1),  and  sQ  *  0. 

For  H.,  we  see  that  (A. 2)  Is  satisfied  with  e  «  0  and  thus  we  impose 
the  restriction  \  <  1.  Then  from  (3.1)  and  (4.1)  we  find 

(4.2)  eu  (Jc,  J)  =  5(3  -  X)3/(27(X2  -  2X  +  5)}. 

"l 

Values  of  (Jc,  J)  are  given  in  Table  4.1.  From  (4.2)  we  note  that, 
as  is  to  be  expected,  as  A  tend:  to  0  (corresponding  to  the  case  of  no  censor¬ 
ing),  eM  (Jc,  J)  tends  to  1. 

H1 

In  order  to  provide  a  reference  point  to  the  amount  of  censoring,  and 

thereby  facilitate  the  interpretation  of  eu  (Jc,  J),  we  also  include  in 

1  -1 

Table  4.1  the  value  of  pH  -  P(X1  <  Yj)  ■  (1  +  A)  ,  the  probability  of 
obtaining  an  uncensored  observation  when  Xj,  Is  exponential  with  scale  para¬ 
meter  1  and  Yj  Is  independent  of  X^  and  has  the  censoring  distribution  Hj. 
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When  the  censoring  distribution  is  Hg,  straightforward  but  tedious 
calculations  yield 


{4’3)  eH2(Jc,  J)  =  ^27(aJ  -  2Ar  +  5)(5crr1(3  -  Arr3[exp{-(3  -  Xr)sf_1> 

m+i  ,  -> 

-  exp {-(3  -  Ar)sr}]  -  S  108(1  -  Ar)(5cr)  *(3  -  \r)  <;[sr_1exp{-(3  -  Ar)sr-1> 

r=l 

m+i  ,  ,  ? 

-  spexp{-  (3  -  Ar)sr}]  +  £  108(5cr)-1(3  -  Arr1[s^1exp{-(3  -  Ar)sf} 

r=l 

-  sJexp{-(3  -  Ar)sr}], 

where  s^  *  •.  Again,  with  the  censoring  distribution  Hg,  (A.2) 

can  Include  the  case  «  ■  0  and  thus  here  we  have  Am  1. 

Values  of  eu  (Jc,  J)  are  also  given  in  Table  4.1,  Again,  as  a  reference 
H2 

point  for  the  amount  of  censoring  under  the  censoring  distribution  Hg,  we 
Include  In  Table  4.1  values  of  pH  =  (Xj  <  Y^)  when  is  exponential  v/ith 
scale  parameter  1,  and  Yj  is  independent  of  X  with  distribution  Hg.  Direct 
calculations  show 
m+1 

pH  ■  1  -  E  crAr(l  +  Ar)  [exp{-(Ar  +  l)sr-1)-  exp{-(Ar  ♦l)sr}],  v/here  sm+1  =  • 
2  r*l 


TABLE  4.1 

Efficiency  loss  under  exponential  (H^)  and  piecewise  exponential  (Hg)  censoring. 

(Hx) 

A:  1  1/2  1/3  1/4  1/10 

eu  (JC,J):  .371  .681  .790  .844  .939 

M1 


500 


667 


750 


800 


.909 
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(h2) 

m=i 

Sl: 

1 

1 

i 

2 

( X  ^ ,  ^  )•’ 

(1/2, 

1) 

(1,  1/2) 

(1/2, 

1/3) 

(1/2,  1) 

eH2^jC»  J): 

.529 

.498 

.723 

.642 

PH  : 

"2 

.630 

.523 

.685 

.675 

m=2 

(sl*  s2^: 

(1/2, 

1) 

(1/2, 

1) 

(1/2, 

1) 

(Xj,  x2,  Xj 

):  (1,  1/2, 

1/3) 

(1/3, 

1/2,  1) 

(1/2, 

1/3, 

1/4) 

.617 

.597 

.772 

V 

.576 

.667 

.718 

5. 

An  Example. 

The  data 

in 

Table 

5.1  are  found  in 

Hollander  and  Proschan 

(1979)  and  are  an  up-dated  version  of  data  given  by  Koziol  and  Green  (1976). 
The  data  correspond  to  211  state  IV  prostate  cancer  patients  treated  with 
estrogen  in  a  Veterans  Administration  Cooperative  Urological  Research  Group 
study.  At  the  March,  1977  closng  date  there  were  90  patients  who  died  of 
prostate  cancer,  105  who  died  of  other  diseases,  and  16  still  alive.  Those 
observations  corresponding  to  deaths  due  to  other  causes  and  those  correspond¬ 
ing  to  the  16  survivors  are  treated  as  censored  observations  (withdrawals). 

As  reported  by  Koziol  and  Green  (1976),  there  is  a  basis  for  suspecting 
that  had  the  patients  not  been  treated  with  estrogen,  their  survival  distri¬ 
bution  for  deaths  from  cancer  of  the  prostate  wuld  be  exponential  with  mean 
100  months. 

Hollander  and  Proschan  (1979)  developed  a  goodness-of-fit  procedure 
for  testing,  in  the  randomly  censored  model,  that  F  is  a  certain  (completely 
specified)  distribution.  They  applied  their  test,  and  competing  procedures 
of  Koziol  and  Green  (1976)  and  Hyde  (1977),  to  the  data  of  Table  5*1*  The 
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hypothesized  F  was  taken  to  be  exponential  with  mean  100.  The  two-sided 
P  values  obtained  were  consistent  with  the  hypothesis.  However,  Gregory 
(1979)  has  proposed  some  goodness-of-fit  tests  which  (for  certain  alterna¬ 
tives)  may  be  more  powerful  than  the  tests  of  Hollander  and  Proschan  (1979), 
Koziol  and  Green  (1976),  and  Hyde  (1977).  Gregory's  tests,  applied  to  the 
data  of  Table  5.1, strongly  indicate  a  deviation  from  the  postulated  expo¬ 
nential,  with  mean  100,  distribution. 

Possible  alternative  models  include  an  exponential  distribution  with 
a  mean  different  than  100,  or  a  distribution,  such  as  an  t©U  distribution, 
that  could  represent  "wearout."  To  explore  the  possiblity  of  the  latter 
type  of  alternative,  it  is  reasonable  to  apply  the  test  based  on  j£. 

Applying  our  NBU  test  to  the  data  of  Table  5.1,  we  obtain 

J211  *  *193*  5211  *  *105  and  (2n)1/2{Jgn  -  (lM)}^  *  -2.56  with  a  cor¬ 
responding  one-sided  P  value  of  .0052.  Thus  the  test  indicates  strong  evi¬ 
dence  of  wearout  and  suggests  that  an  NBU  model  is  preferable  to  an  expo¬ 
nential  model . 


TABLE  5.1 

Survival  times  and  withdrawal  times  in  months  for  211  patients 
(with  nunber  of  ties  given  in  parentheses) 

Survival  times:  0(3),  2,  3,  4,  6,  7(2),  8,  9(2),  11(3),  12(3),  15(2),  16(3), 
17(2),  18,  19(2),  20,  21,  22(2),  23,  24,  25(2),  26(3),  27(2),  28(2),  29(2), 

30,  31,  32(3),  33(2),  34,  35,  36,  37(2),  38,  40,  41(2),  42(2),  43,  45(3), 

46  ,  47(2),  48(2),  51,  53(2),  54(2),  57  ,  60  ,  61,  62(2),  67  ,  69  ,  87  ,  97(2), 

100,  145,  158. 

Withdrawal  times:  0(6),  1(5),  2(4),  3(3),  4,  6(5),  7(5),  8,  9(2),  10, 

11,  12(3),  13(3),  14(2),  15(2),  16,  17(2),  18(2),  19(3),  21,  23,  25,  27, 

28,  31,  32,  34,  35,  37,  38(4),  39(2),  44(3),  46,  47,  48,  49,  50,  53(2), 

55,  56,  59,  61,  62,  65,  66(2),  72(2),  74,  78,  79,  81,  89,  93,  99,  102,  104(2), 
106,  109,  119(2),  125,  127,  129,  131,  133(2),  135,  136(2),  138,  141,  142, 

143,  144,  148,  160,  164(3). 


-  24  - 


Ac know! erf  cement .  Ue  are  grateful  to  Professor  Jayaraman  Sethuraman  for 
many  helpful  discussions. 
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A  life  distribution  F  is  new  better  than  used  ( NBU )  If  F(x  +  y)  £F(x)F(y)  for 
all  x,  y  >  0  (F*  1  •  F).  Using  a  randomly  censored  sample  of  size  n  from  F,  we  propose  a 
<est  of  Hq:  f  is  exponential,  versus  Hji  F  is  NBU,  but  not  exponential.  Our  test  is  based 
on  the  statistic  -  JJFn(x  +  y)dFn(x)dFn(y),  where  Fn  is  the  product  limit  estimator  of  F, 
Introduced  by  Kaplan  and  Meier  (1958). 

Under  mild  regularity  on  the  amount  of  censoring,  the  asymptotic  normality  of 
suitably  normalized.  Is  established.  Then  using  a  consistent  estimator  of  the  null  standard 
deviation  of  an  asymptotically  distribution- free  test  Is  obtained.  Finally,  using 

tests  for  the  censored  and  uncensored  models  we  develop  a  measure  of  the  efficiency  loss 
due  to  the  presence  of  censoring. 


